Pairing in the spin sector 
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The nanometer-scale gap inhomogeneity revealed by recent STM images of BSCCO surface 
suggests that the "gap coherence length" is of that order. Thus a robust pairing gap can develop 
despite the poorly screened Coulomb interaction. This can be taken as an evidence that the pairing 
in high Tc materials occurs primarily in the neutral (spin) sector. We provide theoretical support 
for this point of view. 
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It was first emphasized to us by Pan and Davis 0] 
that the energy gap extracted from STM measurements 
of BSCCO surface is inhomogeneous at nanometer scale. 
Moreover this is true for systems ranging from under- 
to slightly over-doping. |^|-^ Pan conjectured that such 
inhomogeneity is due to the variation of the carrier den- 
sity in the copper-oxygen plane caused by randomly po- 
sitioned dopant oxygen. 

The above findings suggest that the "gap coherence 
length" of high-Tc materials is at most a few nanometer. 
Since at such short length scale the Coulomb interaction 
is poorly screened, it must be true that the pairing takes 
place primarily in the neutral (spin) sector and hence 
hardly affects the charge density correlation. 

Ever since the BCS theory, superconductivity has been 
attributed to the pairing of electrons. In the case of high- 
Tc superconductors, it is sometimes stated that supercon- 
ductivity requires the binding of doped holes. Concep- 
tually it is important to distinguish binding from pair- 
ing. The former is a feature in density-density correlation 
while the latter is manifested in two-particle off-diagonal 
correlation. 

The distinction between pairing and binding is partic- 
ularly pertinent in the cuprates because of the short gap 
coherence length. (By gap coherence length we mean the 
minimum length required for the pairing gap to be ex- 
hibited.) Based on the STM results we argue that such 
length is of order of nanometers. This in turn suggests 
that the pairing responsible for high- Tc superconductiv- 
ity takes place primarily in the spin degrees of freedom. 

We support our point of view by first demonstrating 
that despite strong Coulomb interaction in the following 
Hamiltonian 



H - -t^(c]„c„+ /i.e.) + J^(S, • S,- -n,n^) 
+ Vc'^—{ni-n){nj-n), (1) 



it is energetically favorable for the following variational 
ansatz [Hj to develop d-wave pairing: 



|*(i?,A))=PG|^m/P,A)). (2) 

Here Pq is the projection operator that removes double 
occupancy and |^m/(_D,A)) is the ground state of the 
following mean-field Hamiltonian 



H„, f{D, A)= '^[Xkcl^Cka + iDkc\^Q ^Cko 

k 

+ (Afc4TC-fci+h.c.)], 



(3) 



where Xk ~ — 2(cos/cj; + cos ky), Dk 



cos ky). 



- 2D{coskx — 
In the rest of 



cosky), and Afc = 2A(cosfca; 
the paper we use t/ J = 3. 

By a straightforward minimization we conclude that 
for a wide doping range it is energetically favorable to 
develop a non-zero A (not D) for Vc as big as 9 J. This re- 
sult implies that a non-zero A hardly perturb the density- 
density correlation, hence cannot cause hole binding. We 
emphasize that the purpose of the present discussion is 
not to argue that Eq. (|^) is the ground state of Eq. (|^, 
rather it is to show that the pairing exhibited by Eq. (S) 
is not hole binding. 

To appreciate the effect of the projection operator 
in Eq. (1) we have also investigated the stability of 
4'm/(0, A)) without the Gutzwiller projection. In a 
nearest-neighbour repulsion model the analytic result 
suggests that once Vnn, the nearest- neighbour interac- 
tion strength, is larger than « 0.5J, pairing is absent! 

Now we provide the details. We minimize E{D, A) = 
{'^'\H\'^)/{^\^) by varying D and A. The results pre- 
sented below are obtained for a system with the average 
number of holes per site (x) equal to 12%. The evaluation 
of E{D, A) is performed by Monte Carlo when necessary. 

In Fig. 1(a) we present the projected results for AE = 
E{0, A) - E{0, 0) versus A in a 10 x 10 lattice with 12 
holes. The blue open circles are for the pure t-J model, 
the black crosses are for t-J model with a nearest neigh- 
bor repulsion Vnn = 3 J, and the red open squares are for 
t-J -I- Coulomb model (Eq. (0)) with Vc = 3 J. For each 
of the three cases a nonzero A develops. 

In Fig. 1(b) we present the corresponding results for 
the nearest-neighbor repulsion model when the projec- 
tion operator is removed from the variational ansatz. As 
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one can see even for Vnn = the optimal value of A is 
suppressed. Moreover for > 0.5J the optimal A is 
zero. 

Since the presence/absence of spontaneous staggered 
current order is a timely issue, we have also studied 
the optimum form of Eq. (|^) allowing both D and A. 
In Fig.l(c) we plot AE = E{D,A) - E{0,0) for the t- 



J+Coulomb model {Vc = 3J) at a; = 12%. It is clear 
from this plot that the minimum corresponds to a non- 
zero A but D — 0. Consequently we conclude that for 
the doping relevant to the present paper d-wave pairing 
is the only order that occurs for the model described by 
Eq. §). ' 




FIG. 1. Results for x = 12%. (a) AE as a function of A at D = from Gutzwiller projection in a 10 x 10 lattice. Open circles: 
pure t-J model; crosses: t-J model with nearest neighbor repulsion Vnn ~ 3J; squares: t-J model with Coulomb interaction 
Vc = 3J. (b) Without the Gutzwiller projection, analytic results for nearest-neighbor model, c) AE{D, A) as a function of D 
and A for the projected wavefunction at Vc = 3J. 



The next question is, when subjected to disorder po- 
tential, whether the wavefunction in Eq. (^) can adjust its 
pairing parameter A "adiabatically" to the local density 
to account for the observed gap inhomogeneity. Unfortu- 
nately with the disorder potential breaking the transla- 
tion symmetry the vast variational freedom renders vari- 
ational Monte-Carlo impossible. What we shall do in the 
rest of the paper is slave-boson mean-field theory which 
takes the projection operator in Eq. (H) into account in an 
average fashion. The hope is that since such mean-field 
theory qualitatively captures the physics of Eq. (||) in the 
absence of disorder it will produce meaningful result in 
its presence. 

The starting point of our mean-field theory is the fol- 
lowing Hamiltonian Q] 



H = -tJ2ibMflhc.+ h.c.) + J^(S. ■ S,- 



(4) 



In the above b and / are the holon and spinon oper- 
ators obeying + bjbi = 1, Ui = 1 - bjbi, and 
Si = (1/2) f-^(Jai3fii3. The disorder potential Ui is given 
by 



E 



imp- 



(5) 



In the above Vimy is the location of the oxygen dopant 
and d, the setback distance, is the vertical separation of 



the oxygen dopant plane from the copper-oxygen plane. 
Q From simple valence counting we expect Nim^ to be 
half the number of holes. The results reported in the rest 
of the paper are obtained on 24 x 24 lattice for doping 
X = 12% using Vd/2 = = t ^ 3J, and d = 2x the 
Cu02 lattice constant. We have checked that the results 
change smoothly with d and x. 

The mean-field trial wavefunction 1-0) is given by \tp) — 
\ipb) ® IV'/) where the bosonic (IV'b)) ^^^d fcrmionic (IV"/)) 
states are 



(6) 



n J 



The mean- field single-particle wavefunctions Unj , Vnj and 
Xj are varied to minimize 

i 

-^iJ2{^\n.,~n\^). (7) 

i 

Lagrange multipliers Xi, and fi are introduced to guar- 
antee that the average occupation obeys the constraints 
locally as well as globally. Assuming that the inhomo- 
geneity in STM images are indeed induced by the dopant 
disorder, then, this model should capture the essential 
features of what is seen on BSCCO surface. 
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The biggest difference between screening in an ordi- 
nary metal and in the cuprates is that the latter is close 
to the Mott- insulating limit. Due to the no-double- 
occupancy constraint the ability for charge to redis- 
tribute is severely hindered. For example while the con- 
straint has no effect on the electron depletion, it docs 
forbid local electron accumulation beyond one electron 
per site. As a consequence, in the distribution of holes 
there can be local peaks of hole density while the op- 
posite, i.e. sharp local depletion of holes, tends not to 
occur. This is indeed what comes out of our calculation. 

In Fig. 2 we show the bare impurity potential Ui (Fig. 
2(a)) and the corresponding hole distribution (Fig. 2(b)). 
By comparing the two figures, one sees that the hole dis- 
tribution does correlate with the bare potential. 



a. Impurity Pot. b. Hole Density 




FIG. 2. A plot of the (a) bare impurity potential Ui and 
(b) the hole distribution (b) in a 24 x 24 lattice with x = 12% 
and d = 2. 



We now discuss the mean-field prediction of the tun- 
neling spectroscopy. The local differential conductance 
measured at a bias V by STM is proportional the elec- 
tron local density of states. In the slave-boson theory if 
one writes bi — \bi\(l)i {(p is a U(l) phase factor) and ig- 
nore the fluctuation in the electron spectral function 
is given by 

Pe4V) = hfpgpAV), (8) 

where Pqp^i(V) is the local spectral function of the quasi- 
partides and is the local hole density. (The quasipar- 
ticle creation operator is given by (j)ifia ) I^i the mean- 
field theory where the holon phase fluctuation is ignored 
one obtains 

PqpAV) - \Um\^KV-En)+Y. \Vr,^\^6iV+Er,). (9) 

n n 

In the above (uni, Vni) is the Bogoliubov-deGennes eigen- 
functions of the spinon self-consistent mean-field Hamil- 
tonian. 

From Eq. ^ it is clear that when integrated over en- 
ergy V , the quasiparticle spectral function obeys a sum 
rule, 

/oo 
Pqp,,{V)dV = 2, (10) 
-oo 

hence is independent of the site index i. This is not true 
for the electron spectral function. Due to the presence of 



the factor the total integrated value of the electron 
spectral function depends on the local hole density and 
hence varies from site to site! Such lack of spectral con- 
servation is a generic property of a Mott insulator [ p^ . 
Equation (^) suggests that by properly dividing out the 
integrated electron density of states (and hence \hi\^) one 
can get the quasiparticle density of states. 

It turns out that in the actual experiment this is cus- 
tomarily done. In an STM experiment it is common to 
have the feedback loop set up so that the total tunneling 
current (i.e. the integral of the differential conductance 
up to a voltage Vmax) is held at a fixed value [^j4| This 
way of calibration precisely divides out in Eq. (||). 
Thus the tunneling spectra presented in Ref. |^ is, in 
our language, the quasiparticle local density of states. 
One can also undo the calibration to restore the |6iP and 
hence obtain the electron local density of states g]. It 
turns out that these two density of states have interesting 
observable differences. 

In Figure 3, we plot the quasiparticle {pqp.iiV)) and 
the electron {pf.^i{y)) local density of states in the bias 
range of — 0.6J < < -1-0. 6J for four different locations 
in Fig. 2. Among the four curves the local hole density 
(|6i|^) varies from 0.055 to 0.174. As one can see, the 
peak-to-peak distance, i.e. the local gap, varies consid- 
erably among the curves. (Hj^,^. 

Let' s now focus on the behavior of pqp ,,(t/),and p,,,{V) 
at small V . While different Pqp^i{V) curves tend to merge 
at small voltage, the Pe,i{V) curves do not. This differ- 
ence is precisely caused by the fact that each curve has 
a different \hi\'^ . It is amazing that this difference is seen 
in the experimental curves by changing the calibration! 




-0.6 -0.3 0.3 .^ -0.6 -0.3 0.3 0.^ 



FIG. 3. Tunneling conductance of the fermions {ps,i(V)) 
and of electrons (pe,i(^)) at several positions in the 24 x 24 
lattice (same as in Fig. 1). The local hole densities corre- 
sponding to each curve is indicated in the figure. The bias 
voltage V is in units of J. 

What is the reason for the universality of the quasipar- 
ticle density of states at low energy? The answer is the 
robustness of the low-energy quasiparticle wavefunctions 
due to the vanishing density of states. In Fig. 4 we show 
for two of the low- lying energy eigenstates (associ- 
ated with two different nodes) from two different disorder 
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realization, (a)-(b) and (c)-(d) respectively. These eigen- 
states show a simple, geometric pattern, insensitive to the 
underlying disorder. The orientation of the wavefunction 
is also consistent with the direction of the wavevector 
of the nodal quasiparticles. Such regularity persists up 
to about 0.15J, the same energy scale below which the 
quasiparticle density of states appears universal. 



E/J=0.090 
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FIG. 4. Low-lying eigenstates, lunip, for two disorder real- 
izations ((a)-(b) and (c)-(d)) for several eigenenergies. 



is likely that this discrepancy is due to the omission of 
quantum fluctuation of holon phase (0^) in our calcula- 
tion. 

It has recently been reported that after thermal an- 
nealing, the inhomogeneities of the BSCCO surface dis- 
appears. Such result raises question as to whether the ob- 
served spectral inhomogeneity is an intrinsic bulk prop- 
erty. The point of view we take in this paper is that even 
if the surface inhomogeneity is not intrinsic it still tells 
us an important information, i.e., the existence of a sit- 
uation where the gap varies on nanometer length scale. 
We argue that such nanometer-scale variation suggests 
that the pairing in cuprates occurs essentially in the spin 
sector. 
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Additional comparison can be made between the ex- 
periment and the mean- field result. In Ref. the au- 
thors plot the integrated local density of states vs. the 
local gap. Within the scattering of the data the result 
follows a monotonic trajectory implying a larger gap for 
a smaller and vice versa. Figure 5 is such a plot from 
our mean-field theory. (In making this plot we included 
the local density of states for all 24 x 24 sites in two disor- 
dered samples.) The result agrees with the experimental 
findings qualitatively. 
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FIG. 5. Local hole density vs. local peak-to-peak conduc- 
tance plot. 

Before closing a caveat is in order. One aspect of our 
mean-field result disagrees with the experimental findings 
- our conductance curves for larger gaps show a taller 
peak while the experimental finding is the reverse. It 
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